Weak decays of ^He 
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We measured the lifetime and the mesonic and nonmesonic decay rates of the ^He hypernucleus. 
The hypernuclei were created using a 750 MeV/c momentum K~ beam on a liquid 4 He target by 
the reaction 4 He(_R'~, 7r~)AHe. The ^He lifetime was directly measured using protons from Ap — > np 
nonmesonic decay (also referred to as proton-stimulated decay) and was found to have a value of 
r = 245 ± 24 ps. The mesonic decay rates were determined from the observed numbers of 7r~'s 
and 7r°'s as Y^-/T to t = 0.270 ± 0.024 and IV /T to t = 0.564 ± 0.036, respectively, and the values 
of the proton- and neutron-stimulated decay rates were extracted as r p /r to t = 0.169 ± 0.019 and 
Yn/Ttot < 0.032 (95% CL), respectively. The effects of final-state interactions and possible three- 
body ANN decay contributions were studied in the context of a simple model of nucleon-stimulated 
decay. Nucleon-nucleon coincidence events were observed and were used in the determination of the 
nonmesonic branching fractions. The implications of the results of this analysis were considered for 
the empirical AI = | rule and the decay rates of the hypernucleus. 

PACS numbers: 21.80.+a,13.75.Ev 



I. INTRODUCTION 

Studies of the hyperon-nucleon interaction investigate 
weak processes that are not well understood in multiple 
baryon systems. In particular, the study of nonmesonic 
decays of light hypernuclei is the only tractable method 
of investigating the AS = 1 weak baryon-baryon inter- 
action. By comprehensive measurements in s-shell hy- 
pernuclei, various A7V initial states can be investigated. 
Few-body nuclear structure calculations can then relate 
the decay rates to the underlying weak AN interactions. 
These studies are complementary to investigations of the 
parity- violating contributions to the nucleon-nucleon in- 
teraction. However, the AN system can be prepared in 
a particle-stable state so that only weak processes con- 
tribute to the decay rate. 

A particularly interesting aspect of AS = 1 transitions 
is the empirical AI = \ rule. We note that if strong inter- 
actions are neglected, the AS = 1 non-leptonic decays of 
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hyperon systems would be expected to be dominated by a 
weak quark-level Hamiltonian based on one W-exchange 
of the form: 

G 

#weak = —j= sin C cos 6 C (u7 P (l - 7 5 )s(V'(l - 75 )w) . 

The isospin character of this fundamental Hamiltonian 
does not favor AI = | weak transitions over the AI = | 
ones. However, it has long been recognized that the 
AI = i transitions dominate all observed non-leptonic 
decays of strange hadrons. Although it is generally be- 
lieved that the presence of the strong interaction accounts 
for this empirical AI = ^ selection rule, the various expla- 
nations for the dominance of the AI = \ amplitudes [TJ |2] 
have indicated that the mechanism may be specific to the 
given system and not a universal feature of the AS* = 1 
weak Hamiltonian. It remains an open question whether 
all AS = 1 non-leptonic decays obey the rule or it is spe- 
cific to the pionic decay channels (the only tested class of 
non-leptonic decays). The study of hypernuclear decays 
provides an opportunity to answer this question by pro- 
viding the only experimentally accessible non-pionic test 
of the AI = 1 rule. 

The total decay width of a hypernucleus includes con- 
tributions from mesonic and nonmesonic decay modes. 
The mesonic decay modes, defined experimentally by the 
observation of the appropriate pion in the final state, 
are dominated by single-body processes analogous to 
the decay of the free A, A — > 7r~p and A — > 7r°n, 
along with additional multi-body contributions such as 
Ap — » nmr + . The nonmesonic modes are expected to be 
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dominated by the proton-stimulated reaction Ap — > rap 
and neutron-stimulated reaction An — > rara with rates T p 
and T„, respectively. Other multi-body interactions such 
as ANN — > NNN may also contribute with rate V m b- 
Determination of the nucleon-stimulated decay rates is 
more complicated than the mesonic case and is ultimately 
model dependent. The total decay width can be written 
in terms of the widths of these mesonic and nonmesonic 
decay modes as: 

f\otal r mes0 nic r nonmGSO nici 

. > , ' > (2) 

— iv- + f ot d + rjr+ + r p + r„ + r m &, 

where contributions from semileptonic and weak radia- 
tive A decays have been neglected as they contribute only 
about 0.3% to the total free A decay width jjJJ. 

The two-body nonmesonic decay modes, AA" — > 
NN, are readily distinguishable from the mesonic decay 
modes, A — ► nN, because of the large energy (M\—M n — 
176 MeV) available for the final-state nuclcons. These 
modes are sensitive to weak interaction couplings (such 
as g\Np or gNNi<) not available to the free hyperon de- 
cays. Also, it has been suggested that several of the 
weak AN — > NN amplitudes are dominated by direct- 
quark processes 4 5, 6 in which no intermediate meson 
is present in the interaction. Such direct processes may 
not adhere to the A I — ^ rule. 

In nuclei the A-particle is not Pauli blocked, thus hy- 
pernuclei generally decay with the A-particle initially in 
the lowest ls-state. In order to determine properties of 
the underlying AA~ — > NN interaction, it is convenient 
to consider systems that confine the nucleons to the low- 
est orbitals. Figure [TJ shows the spin-isospin character of 
the initial (AN) and final (nN) states for s-shell hyper- 
nuclei in a single-particle shell model picture. The 1=1 
final states are accessible to both proton- and neutron- 
stimulated decays, while the 1=0 final states are avail- 
able only for proton-stimulated decay. Further complica- 
tions arise as one increases A because the high probability 
of final-state interactions dilutes the correlation between 
nucleons observed and nucleons responsible for the decay. 
Table [I] shows that a complete set of measurements of the 
nonmesonic decay widths for the s-shell hypernuclei al- 
lows one to isolate specific AN initial states. Thus, one 
can determine the spin-isospin structure of the funda- 
mental AA" — -> nN weak interaction by combining mea- 
surements of these few-body hypernuclei with detailed 
finite nucleus calculations that take into account the dif- 
ferences in the initial and final state phase-space. How- 
ever, at this time the experimental data are incomplete 
and often have large estimated errors. 

A phenomenological model first put forth by Block and 
Dalitz [5] can be used to relate the nonmesonic decay 
rates of the s-shell hypernuclei listed in Table [TJ Specif- 
ically, the current experimental results for ^He (this 
work) and iHe (Outa et al [9]) can be extended to the 
iH hypernucleus, for which measurements are currently 
scarce. The model employs the mean nucleon density, 



TABLE I: Allowed initial states for Is— shell hypernuclei. 
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FIG. 1: Spin-isospin selection rules for the AN — > nN process 
in ls-shell hypernuclei. 

pA, at the location of the A and makes several assump- 
tions concerning the AA~ interaction: (1) the nucleon- 
stimulated A decay is treated as incoherent, (2) final- 
state interactions are not included, and (3) three-body 
nucleon-stimulated decays are neglected. These assump- 
tions are seen to be adequate at the present level of accu- 
racy of the experimental quantities. Defining the quan- 
tity Rns as the rate for an A^-stimulated decay originat- 
ing from AN relative spin state S, the nonmesonic decay 
rates for the s-shell hypernuclei iH, iHe, and iHe are 
derived as: 

r„ m (iH) = 7T (RnO + 3i?nl + 2-Rpo) j (3) 

6 

r nm (iHe) = -— - (2-Rno + Rpo + 3.Rpi) , (4) 
r nm (^He) = ~ (R n o + 3i?„i + R p o + 3R p i) , (5) 

o 

where the mean nucleon density, pA, is defined as: 

PA = J PA(r)\^(r)\ 2 dr, (6) 

for the nucleon density, pa(^), and the A wave function, 
ip(r). Also, it has been assumed that the mean nucleon 
density, P4, has the same value for both iH and iHe. 

Taking the ratio of the rates of neutron-stimulated 
iHe decay to proton-stimulated iH decay yields: 

r n (jHc) = Rno = 2 ^ 
Tp(aH) R P o 

where the value of 2 is the Al = | rule prediction for 
this ratio. Similarly, an expression for the iH neutron- 
stimulated rate can be found by considering the ra- 
tio of neutron-stimulated iH decay to proton-stimulated 
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^He decay: 



r n (iH) _ Rr, 
r p (4 He ) ~ r 



pO 



3R n 
~3R 



r„ 



(iHe), 



(8) 



independent of any assumptions about the AI = g rule. 

In this article, we present the results of BNL exper- 
iment E788, which measured the lifetime and partial 
widths of the \He hypernucleus. Section [TT] presents the 
experimental apparatus used for these measurements, 
and a discussion of the out-of-beam tracking and par- 



ticle identification algorithm is given in Section III The 
excitation energy spectra for the 4 He(i"C~, 7r~)^He pro- 
duction reaction are discussed in Section HVl The extrac- 
tion of the lifetime of \Rc and particle-emission spectra 
are presented in Sections [V] and VI respectively. The 
determination of the ^He mesonic and nonmesonic par- 
tial decay rates is discussed in Section |VII| In the case 
of the nonmesonic proton- and neutron-stimulated de- 
cays, the rates are found from both the single-particle 
kinetic energy spectra of Section |VI| as well as from ob- 
served multiple-nucleon coincidence events, giving two 
essentially independent determinations. A comparison 
with the results of other experiments and theoretical 
calculations and the extension of these results to the 



i H hypernucleus is presented in Section VIII 



II. EXPERIMENTAL SETUP 

The partial decay rates and lifetime of ^He were mea- 
sured using the LESB-II kaon beamline of the AGS at 
Brookhaven National Laboratory. The beamline param- 
eters are given in Table [IT] The experimental layout is 
shown in Fig. [2] A beam of 750 MeV/c kaons inci- 
dent on a liquid helium target produced ^He through 
the 4 He(X~, ir~)\He reaction, and pions from the pro- 
duction reaction were detected near 0° in the Moby Dick 
spectrometer. Detector packages located on each side 
of the target area were used to detect and identify pions, 
protons, neutrons, and 7's from the hypernuclear decays. 

The incident kaon trajectories were determined with 
drift chambers located just upstream of the target (ID2X 
through ID7U). This information was combined with 
data from the drift chamber ID IX, located upstream of 
the K~ spectrometer, to determine the incident kaon 
momentum. Although the kaon-to-pion ratio at the 
target was only 1:13, the kaons were cleanly sepa- 
rated using a combination of the time-of-flight between 
scintillator IS1 (not shown) located upstream of the 
K~ spectrometer and scintillators IT1 and IT2 and a 
critical angle Cerenkov detector (CP). The outgoing pion 
momentum was measured with the Moby Dick spectrom- 
eter located at zero degrees downstream of the target, 
and the outgoing particle identification was accomplished 
with time-of-flight measurements. 

The 30.5 cm long liquid helium target vessel con- 
sisted of a 6.35 cm diameter Dacron cylinder with a wall 




Beam Left 



K- Spectrometer 



Beam Right 



FIG. 2: Shown here is the E788 experimental apparatus 
used to carry out the ^He hypernuclear decay measurement. 
Also shown are the focusing quadrapoles at the exit of the 
K~ spectrometer and at the entrance of the Moby Dick spec- 
trometer. 



TABLE II: LESB-II and Moby Dick Spectrometer. 

Item Description 

LESB-II Separated Beamline/Spectrometer 
Momentum 750 MeV/c 

Momentum acceptance 5% FWHM 

7r : K ratio at target 13 : 1 

K~ flux 2 x lOV 1 

Momentum measurement Drift chambers ID1X - ID7U 

and hodoscope MH 
7T - K Particle ID IS1 - IT1.2 TOF scintillators 

and CP Cerenkov 7r-veto 



Moby Dick Spectrometer 



Momentum 

Momentum acceptance 
Momentum measurement 
7r — K particle ID 



680 MeV/c 

5% FWHM 
Drift chambers ID9Y - ID15Y 
IT1,2 - IS2,3 TOF scintillators 



thickness of 0.020 cm. To minimize material between 
the decay vertex and the out-of-beam decay particle de- 
tectors, the vacuum vessel containing the Dacron shell 
was constructed from a PVC foam cylinder of density 
0.053 gm/cm 3 with inner and outer radii of 7.6 cm and 
11.43 cm, respectively. 

The separation of the decay protons and pions was 
accomplished by combining the crude range information 
from all out-of-beam detectors with the measured rate 
of energy loss (dE/dx), total energy loss, and the mea- 
sured time-of-flight (TOF) between the Decay-Timing 
and TOF scintillator layers (also referred to in this pa- 
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per as 0T1 and OT2, respectively). The precision tim- 
ing scintillators of the Decay-Timing layer also played 
an important part in the determination of the lifetime 
of the ^He hypernuclcus. With no magnetic field in 
the region of the out-of-beam detectors, the detector 
system of E788 was unable to distinguish a n~ from a 
7r + , making it possible to measure only the sum of the 
two widths r„— and T 7r +. However, the width 1"^+ has 
been seen to be small compared to the ir~ decay width 
(IV /IV = 0.043 ± 0.017 [10]) and was ignored in this 
analysis. Large volume (« 1.4 m 3 ) neutron detector ar- 
rays were placed near the target to directly measure the 
decay neutrons and were used to determine the neutron 
energies by time-of- flight. The neutron detector arrays 
were also used to detect the 7's from the decay of 7r°'s 
originating from ^He decay. Thus, this experiment di- 
rectly measured r to tai, r w -, T^o, T p , and r„, along with 
the kinetic energy spectra for the decay 7r _ 's, protons, 
and neutrons. Details of the detector dimensions are 



given in Table III 



III. OUT-OF-BEAM TRACK 
RECONSTRUCTION AND IDENTIFICATION 

The task of the out-of-beam tracking was to sort 
through the scintillator hits recorded in the out-of-beam 
detector arrays and determine the trajectories of the 
charged and neutral particles that created them. A hit in 
the out-of-beam detector array was defined as a coinci- 
dence of signals in the photomultiplier tubes at each end 
of a scintillator element, requiring ADC and TDC infor- 
mation for both of the PMTs. Each hit was characterized 
by a hit time, energy deposition, and hit position. For a 
hit to be considered for tracking, the recorded time (mea- 
sured relative to the event start-time as set by the mean 
time of in-beam scintillators IT1 and IT2) and measured 
energy deposition must have been greater than zero. 

As a preliminary step, neutron detector hits that sat- 
isfied the rudimentary conditions described above were 
examined, and correlated hits were grouped into clusters 
based on relative hit times and positions (i.e., groups of 
hits that were sufficiently close in time and space were 
considered to have arisen from a single particle) . The hit 
clusters were classified as charged if the group included 
any hits in the OT2 or NV layers, otherwise the cluster 
was considered to have originated from a neutral particle. 

Charged tracks were indicated by a coincidence of hits 
in the OT1 and OT2 layers. The charged tracking began 
with a search for hits in the OT1 layer, and for each OT1 
hit found, a search was made for corresponding hits in 
OT2. If the hit pairs were properly time ordered (i.e., 
the OT1 hit preceded the OT2 hit), they were considered 
to be part of a charged track and were then subjected to 
a more detailed tracking procedure. The final track was 
found from a linear least squares fit to the OT1 and OT2 
hits, the Kn reaction vertex, out-of-beam drift chamber 
hits, and any hits in the neutron veto layer or neutron 



detector array that belonged to the same cluster as the 
OT2 hit. Interesting quantities such as the particle's ve- 
locity between the OT1 and OT2 layers and the particle's 
range and total energy loss in the out-of-beam detector 
system were then determined. 

The tracking information found above was used to dis- 
tinguish protons from pions, a task that was complicated 
by the absence of a momentum-analyzing magnetic field. 
To determine the identity of the charged particles, the 
measured velocity, range, and energy loss of a charged 
track were compared with the expected behavior of the 
hypernuclear decay products via the quantity: 



^AE meas — AEx^piPfit) 



ftmeas ftfit 
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, (9) 



where AE meas is the measured energy deposition in the 
OT2 scintillator layer, f3 meas is the velocity fraction as 
measured in the time-of-flight region of the out-of-beam 
detector array, and oe ~ l/v£ and <jp ~ 1 are the res- 
olutions for the measured quantities. To find the quantity 
(3 fa (and thus, the theoretical energy deposition AE n p ), 
the above expression for -D^p was minimized by varying 
the value of (3 fa within limits determined by the mea- 
sured range of the particle. The PID (particle identifica- 
tion) for a given track was then determined by testing the 
values of (calculated assuming the track was created 
by a pion) and D p (calculated assuming the track was 
created by a proton) against predetermined limits. The 
particle was then classified as either a proton, a pion, or 
an unknown (neither nor D p was within the appro- 
priate limits). The unknown particles consisted mostly 
of electrons and were well separated from the pions and 
protons. 

A neutral particle appeared as a single hit or small clus- 
ter of hits in the out-of-beam detector arrays, providing 
little information for tracking. A neutral track was taken 
as a straight line connecting the Kit reaction vertex with 
the out-of-beam neutral cluster, where the time and posi- 
tion of the reaction vertex were known from the in-beam 
tracking. For each of the neutral clusters found in the 
preliminary grouping, a cluster hit time, position, and 
energy were calculated. The hit time and position were 
determined by the member hit that was closest to the tar- 
get, and the energy was taken as the sum of the energy 
deposition for all member hits. Additionally, each cluster 
was characterized by a cluster size defined simply as the 
total number of member hits. Once the neutral track was 
determined by the production vertex and cluster hit po- 
sition, the velocity fraction of the neutral particle could 
be calculated. 

The next task was to separate neutrons from 7's and 
attempt to suppress the large accidental neutral back- 
ground. To reduce accidentals, an energy deposition 
threshold of Ed ep > 5 MeVee (MeV electron equivalent) 
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TABLE III: Descriptions of the target area detectors utilized for the hypernuclear decay measurement. An overview of the 
layout of these detector elements is shown in Fig. [2] 

Name Function Dimensions (cm ) 

Drift Chambers 

ID K - 7T trajectories 12.2 x 12.2 x 0.43 

OD y and z track positions (1 pair each side) 30.5 x 30.5 x 1.9 

Scintillation Counters 

IT1, IT2 K - 7T TOF and ^He formation time 15.0x4.0x1.27 

CP critical angle pion Cerenkov counter 15.0 x 4.0 x 0.6 

OT1 decay time and dE/dx (5 each side) 22.0 x 4.5 x 0.64 

OT2 TOF, dE/dx, and range (10 each side) 150.0 x 12.0 x 1.5 

NV range and charged particle veto (9 each side) 182.9 x 22.9 x 0.95 

NA charged particle veto (2 each side) 182.9 x 22.9 x 0.95 

ND range and neutral particle detection (50 each side) 182.9 x 15.2 x 5.1 



was applied to the neutral clusters, and clusters with cor- 
responding hits in the auxiliary veto counters NA were 
discarded. Good neutron candidate events were then se- 
lected by requiring the neutral cluster to have fewer than 
3 member hits and a value for the neutron kinetic en- 
ergy (as derived from the measured time-of-flight) from 
50 to 200 MeV. The lower bound at 50 MeV avoids the 
energy region where accidentals dominate, and the upper 
bound at 200 MeV, which is well above the maximum ex- 
pected kinetic energy for neutrons from hypernuclear de- 
cay (~170 MeV), removes 7's from consideration. Good 
7 events were selected by requiring the measured f3 of 
the candidate track to be in the range 0.666 < f3 < 2.0, 
where the lower bound at 0.666 is well above the highest 
j3 expected for neutrons from hypernuclear decay. Even 
with the cuts described above, a sizeable background still 
remained in the neutron and 7 samples. The subtraction 
of this background will be discussed in subsequent sec- 
tions. A cancellation effect inherent in the subtraction of 
the background reduces the sensitivity of the final results 
to the selection of the cut parameters given above. 



IV. iHe EXCITATION SPECTRA 

Reconstruction of the K~ and ir~ trajectories through 
their respective spectrometers allowed the calculation 
of the invariant mass of the unobserved strange sys- 
tem. Unlike p-shell hypernuclei, the ground state is also 
a substitutional state and thus the [K~ , it~ ) reaction 
near zero degrees is ideal for creating ^He. The main 
source of background for the production reaction was in- 
flight kaon-decays. Those decays that may satisfy the 
(K~,ir~) trigger condition include (i) K~ — > tt~it , (ii) 
K~ — > fM^Pfj,, and (iii) K~ — ► /i~ £^7r° . For kaons that de- 
cayed within the area of the target, the two-body decays 
were easily removed using the decay kinematics, while 
the three-body decays had to be modeled with a Monte 



Carlo simulation and fit to the data. Some kaons de- 
cayed within the bcamlinc chambers or Moby Dick spec- 
trometer or traversed the experimental area completely. 
These events were characterized by a small apparent 
Ktt scattering angle and were removed by requiring an 
angle greater than 25 mrad. Events with production ver- 
tices outside the target region or poor beam track recon- 
struction were also rejected. 

The resulting excitation energy spectrum for unbiased 
(K~,ir~) events is shown in Fig. |3ja), where unbiased 
refers to the trigger condition in which only an incoming 
K~ and outgoing 7r~ were required (i.e., no information 
from the out-of-beam detectors was considered). Also 
shown are the excitation energy spectra for events with 
a coincident proton, pion, or neutron [Figs. [3jt>) , (c), 
and (d), respectively] where the coincident decay par- 
ticles were identified as described in the previous sec- 
tion. (Note that the low number of counts in the un- 
biased spectrum as compared to the coincident spec- 
tra is a result of a hardware prescale in the unbiased 
trigger.) Since the coincident protons were required to 
have kinetic energies greater than 35 MeV, only events 
in which the A was bound to a nucleus and decayed 
through a nonmesonic channel contribute to the proton- 
tagged spectrum of Fig. 3pb). A fit of a Gaussian shape 
to this peak determined an experimental resolution of 
2.3 MeV RMS and defined the zero of excitation en- 
ergy. The unbiased spectrum, Fig. |3ja), was then de- 
termined to consist of three parts: (i) a broad back- 
ground from K~ — > /i~i7(7r contamination whose shape 
was determined by Monte Carlo simulation and whose 
magnitude was found from a fit to the lower part of the 
spectrum, (ii) events corresponding to the production of 
the bound ^He ground state with mean and width deter- 



mined from Fig. |3Jb) and amplitude determined by a fit 
to the left half of the peak, and (iii) additional structure 
above zero excitation energy that we associate with the 
production of unbound A hyperons. 
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Excitation Energy (MeV) 

FIG. 3: The excitation energy spectrum for 4 He(K~ , 7r~)X 
is shown for several cases: (a) unbiased by out-of-beam tags, 

(b) tagged by a proton in the out-of-beam detector package, 

(c) tagged by a pion, and (d) tagged by a neutron. Overlaid 
on each plot is the result of the fit to the ground-state peak 
(including appropriate backgrounds as described in the text). 



tex (and decay point) was inferred from the in-beam 
and out-of-beam tracks using a distance of closest ap- 
proach (DCA) algorithm. The velocity of the incoming 
K~ was derived from its momentum as measured by the 
upstream K~ spectrometer, and its trajectory was deter- 
mined by the six drift chambers located just upstream of 
the 4 He target. A time measurement was also made for 
the beam kaon as it passed the in-beam scintillators IT1 
and IT2. This measurement set the start time for the 
event and was taken as the average of the two counters 
(producing a combined resolution of ctit = 44.3 ps). The 
time of hypcrnuclcar formation was then simply found 
from the velocity of the K~ and its path length from the 
point midway between the IT1 and IT2 timing counters 
to the reaction point. (The velocity of the beam kaon 
was assumed to have remained constant as it traversed 
the experimental area.) For the out-of-beam charged 
track, the energy was determined from the time-of-flight 
as measured between the Decay-Timing (OT1) and TOF 
(OT2) scintillator layers. The decay time at the vertex 
position was determined after making corrections for the 
unmeasured energy losses in the target region. 

The events used in the \Rc mean lifetime determina- 
tion were subject to the in-beam tracking cuts described 
in Sec. |IV| for the hypernuclear excitation energy spec- 
tra. Protons were chosen in the out-of-beam to avoid 
the contamination of quasifree A decay events present in 
the 7r~-tagged sample. Also, protons with kinetic energy 
below 70 MeV were not considered for the lifetime mea- 
surement since the uncertainty in the energy loss cor- 



The excitation energy spectrum for coincident de- 
cay pions of Fig. 3jc) is shown with a two parame- 
ter fit to the left half of the ground-state peak. The 
fit determines the amplitude of the ground-state peak 
and includes a flat background term. The requirement 
of a decay-pion tag virtually eliminates the residual 
K~ — > /i^D^tt contamination but retains events from 
both bound and unbound A production as would be ex- 
pected. For the neutron-tagged excitation energy spec- 
trum of Fig. |3jd), the coincident neutrons consist of a 
mixture of neutrons from nonmesonic ^He decay, acci- 
dental neutrons, and neutrons created when 7r - 's from 
^He decay or quasifree A decay were absorbed in the ma- 
terial around the target area. 



The lifetime of a ^He hypernucleus was defined event- 
by-event as the time between the formation of the hy- 
pernucleus and its subsequent decay as determined from 
the velocity and trajectory of the incident kaon and the 
energy and trajectory of the charged hypernuclear de- 
cay product. The ^He hypernucleus was assumed to be 
created at rest, and the position of the reaction ver- 




0.5 
Time (ns) 

FIG. 4: Figure (a) shows the ^He lifetime distribution for 
events with a coincident proton overlaid with the fit of the 
expected distribution. Figure (b) shows the time-difference 
distribution for protons from the reaction 4 He (tt~ ,p)X used 
to determine the time resolution function. 
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rection increases at lower energies. To ensure that the 
kaon and proton tracks were correlated with a hyper- 
nuclear formation event: (i) the excitation energy was 
required to be within ±2cr of the mean, (ii) the K~ and 
proton tracks must have had a DCA of less than 2 cm, 
and (iii) the reaction vertex must have been within the 
liquid 4 He containment vessel. The resulting lifetime dis- 
tribution for these proton-tagged events can be seen in 
Fig.ga). 

The mean lifetime of ^He was determined by a fit to 
the distribution of Fig.[4^a). The actual time distribution 
expected for this measurement has the form: 



L T (t) 



dt'R(t')P T (t-t'), 



where P T (t) is the decay probability distribution: 



Pr(t) = -e 

T 



(10) 



(11) 



and R(t) is the resolution function as determined from 
the time distribution of prompt protons originating from 
the reaction 7r~ + 4 He — > p + X [shown in Fig. Qb)] and 
9{t) is the usual theta function. The measured distribu- 
tion ofFig.gb) gives a resolution of a r ~ 150 ps. The fit 
to the hypernuclear lifetime distribution of Fig. [4^a) was 
performed using a likelihood function derived assuming 
Poisson statistics from the function: 



f(t) = Cl + C 2 t + ClL T (t), 



(12) 



which includes two background terms in addition to the 
expected lifetime distribution of Eq. |10| By maximizing 
the likelihood function on the four parameter space (C\, 
C-2, Cs, r), the value for the mean lifetime of ^He was 
found to be r = 245 ± 24 ps. [The values for the back- 
ground terms C\ and C'2 were found to be negligible com- 
pared to Cf L T (t).] 



VI. PARTICLE EMISSION SPECTRA 

The kinetic energy spectra for the decay protons, pi- 
ons, and neutrons are shown in Figures [5j a), |5jb), and 
[6jb) , respectively. These events were subject to the cuts 



described in Sec. 



IV 



for the \\ 



v He excitation energy spec- 
tra. To reduce contamination from out-of-beam particles 
not originating from a hypernuclear ground-state event, 
a cut was also placed on the value of the excitation en- 
ergy. In the case of proton- and neutron-tagged events, 
the excitation energy, EE, was required to be within ±2<r 
of the mean of the ground-state peak, while in the case 
of pion-tagged events, a tighter cut of /i — 2er < EE < fi 
(referred to here as the "— 2er cut" ) was employed to re- 
duce contamination from quasifree A decay events. (The 
quantities p, and a are the previously determined mean 
and width of the ^He ground-state peak.) For the proton 
and pion spectra, corrections for the unmeasured energy 
losses in the target region have been included. 
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FIG. 5: Shown here are the kinetic energy spectra for 
(a) protons and (b) pions detected in coincidence with 
^He formation. These spectra have been corrected for PID 
inefficiencies, detector acceptance, and unmeasured energy 
losses in the target region. 



In order to obtain the neutron kinetic energy spectrum 
of Fig. [6jb) , a considerable background consisting of acci- 
dental neutrons and neutrons from tt~ interactions had 
to be subtracted from the raw neutron kinetic energy 
spectrum shown in Fig. |6ja). The expected shape of 
the accidental neutron spectrum was generated assum- 
ing a flat distribution in f3~ l (i.e., a flat time distribu- 
tion). The distribution for neutrons from 7r~ absorption 
was then extracted from the data by considering neutron- 
tagged events with excitation energy in the range 10 < 
EE < 30 MeV. The neutrons from these events consist 
of a mixture of accidental neutrons and neutrons from 
7r~ absorption only. The normalizations for the back- 
ground shapes were found via a detailed comparison of 
the unbiased and neutron-tagged excitation spectra of 
Figures |3ja) and |3jd). The background contributions 
thus determined are overlaid on the plot of Fig. |6ja). 

The kinetic energy spectra of Figures |5ja) , |5jb) , and 
|6jb) also include corrections for small PID inefficiencies 
and detector acceptance which were derived with the aid 
of a Monte Carlo simulation of the E788 out-of-beam de- 
tector system. Protons, neutrons, and 7r~'s with realistic 
initial energy distributions were propagated through the 
target area and out-of-beam detector system. The re- 
sponse of the scintillator counters was modeled using an 
approximate form of the Bethe-Bloch equation of energy 
loss in the case of charged particles and the DEMONS 
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FIG. 6: Figure (a) shows the raw neutron kinetic energy spec- 
trum with background contributions from accidental neutrons 
(dashed line) and neutrons from n~ absorption (dotted line) 
overlaid. Figure (b) shows the neutron spectrum after back- 
ground subtraction and corrections for detector acceptance. 



software package [IT] for neutron interactions. Unmea- 
sured energy losses in the target materials and detector 
wrappings, pion decays and interactions, and corrections 
for the non-linear response of the scintillator material at 
large energy deposition were also included. The tracking 
and particle identification was performed with the same 
code as used for the real data, and the results were used 
to derive the necessary corrections. 



VII. EXTRACTION OF PARTIAL RATES 

The decay modes considered in this analysis were 
the 7r°, 7r~ , proton-stimulated, neutron-stimulated, and 
three-body nuclcon-stimulated modes with branching 
fractions defined as: 



B a = 



r„ 



JV HN 



(13) 



where a = tt° ,tt~ ,p,n, and mb\ N a is the total number 
of hypernuclear events decaying to mode a; and N^ = 
27800 ± 1500 is the total number of ^He hypernuclear 
events found from the fit to the ground-state peak of the 
unbiased excitation spectrum of Fig.[3ja). This total also 
includes corrections for a hardware prescale in the unbi- 
ased trigger. The 7r° branching fraction was extracted 
by considering 7-tagged hypernuclear formation events, 
while the 7r~ branching fraction was determined from 



the observed ir~ kinetic energy distribution. The non- 
mesonic decay branching fractions were extracted by two 
essentially independent methods utilizing (i) the single- 
particle kinetic energy spectra with final-state interac- 
tions and three-body ANN decay contributions consid- 
ered within a simple model, and (ii) multiple nucleon co- 
incidence events. The final values for the ^He branching 
fractions were then determined by a \ 2 minimization that 
considers the mesonic and nonmesonic rates simultane- 
ously. 



A. The tt° Decay Mode 

The 7r° decay rate was found by considering events for 
which one or more 7's were detected in coincidence with 
^Hc hypernuclear formation. The 7's were identified as 



described in Sec. Ill and the candidate events were sub- 
ject to the same cuts as the excitation spectra of Sec. |IV| 
Three possible sources of 7-rays were considered for the 
7r° rate extraction: 7r° decay, 7r _ interactions, and ac- 
cidental hits. Thus, the possible backgrounds include 
7-tagged events arising from accidentals, 7r~'s and 7r°'s 
from quasifree A decay, 7r~'s from ^He decay, and 7r°'s 
from the K~ — > [i~v il Tt° decay background. To find the 
7r° branching fraction, these background events were sub- 
tracted as described below. 

In order to determine the quantitative contributions 
of each type of background, the various sources of 7's 
were isolated by considering 7-tagged events with par- 
ticular excitation energies. First, events with excitation 
energy below —10 MeV, which consist almost exclusively 
of accidental 7's and 7's from the K~ decay background, 
were used to estimate the detection efficiencies for ac- 
cidentals and 7's from ir° decay. Next, the 7's from 
7r~ interactions were studied using events with excita- 
tion energy in the range 10 < EE < 30 MeV. These 
events consist of accidental 7's and 7's from the K~ 
decay background (as in the previous case) with addi- 
tional background 7 events originating from quasifree A 
decay. Using the accidental and ir° decay 7 detection 
probabilities determined above, the contribution from 
quasifree A — * p7r~ decay could be extracted, allowing the 
determination of the detection probability for 7's from 
7r~ interactions. Finally, the y\He tt° decay events were 
selected by requiring the excitation energy to be within 
±2cT of the ground-state peak and subtracting the back- 
ground contributions found above. 

After properly accounting for the backgrounds and 7 
detection efficiency, the ir° branching fraction was found 
to be B^o = 0.552 ± 0.076(stat.) ±0.061(syst.). So, for a 
A embedded within the ^He hypernucleus, the n° decay 
rate: 

^ = B^o x = 0.59 ±0.10, (14) 

expressed here in units of the total decay rate of the free 
A, is enhanced by a factor of about 1.6 as compared to the 
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free A decay partial width of V{ r ee /V K = 0.358±0.005 |3J. 
(The quantity T to t/^A = 1.07±0.11 was derived from the 
A He hypernuclear lifetime found in Sec. |v|. 



B. The 7r Decay Mode 

The 7T _ decay rate was determined by a fit to the 
observed 7r~ kinetic energy spectrum using a theoret- 
ical distribution provided by the model calculation of 
Kumagai-Fuse, et al. |12j Their predicted tt~ spectrum 
[shown in Fig. ma)] was derived using the resonating 
group method and includes final-state interactions and 
pion distortion effects. The effects of detector acceptance 
and resolution unique to E788 were applied to the model 
by using the spectrum of Fig. |7ja) as input to the Monte 
Carlo simulation described in Sec. |VI| and reconstructing 
the initial kinetic energy distribution with the analysis 
code. 
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FIG. 7: Shown here are (a) the theoretical tt~ kinetic energy 
spectrum of Kumagai-Fuse, et al. [TJ] and (b) a fit of this ex- 
pected shape to the measured distribution after including the 
effects of detector acceptance and resolution. The measured 
spectrum of figure (b) has been corrected for PID inefficiencies 
and unmeasured energy losses (but not detector acceptance). 

A fit of the model prediction to the PID-corrected 
7r~ spectrum of Fig. |7jb) was then performed using the 
fit function: 



N^E) = Mfh w (E), 



(15) 



where h„(E) is a parameterization of the expected 
7T _ kinetic energy distribution as reconstructed by the 



analysis code (normalized to the detector acceptance and 
in units of counts/bin per event) and M^ a is the fit pa- 
rameter representing the acceptance-corrected number of 
pions within the —la excitation energy cut. The result 
of the fit is shown in Fig. [7](b) where the value of M% a 
was found to be 3690 ±190. Correcting for the excitation 
energy cut yields a value for the total number of 7r~'s: 



M. 



total 



9v Mf = 7490 ± 480, 



(16) 



where — 2.030 ± 0.076 is the ratio of the number of 
counts in the ground-state peak of the 7r _ -tagged excita- 
tion spectrum of Fig. |3^c) to the number of events within 
the — 2a cut. 

The 7r _ -decay branching fraction for the ^He hyper- 



nucleus was then found by Eq 
0.022(stat.) ± 0.014(syst.) 
decay, the observed ir~ decay rate: 



13 as B„- = 0.269 ± 



In the present case of ^He 



r A 



B„- x 



r, 



0.289 ± 0.037, 



(17) 



in units of the free A decay, is seen to be suppressed 
by a factor of about 2.2 as compared to the free A decay 
value of r^ ee /r A = 0.639±0.005 [3]. This suppression is 
mainly due to the Pauli blocking of the final-state proton. 



C. Nonmesonic rates from single-particle KE 
spectra 

The rates for proton- and neutron-stimulated decay 
were extracted simultaneously from the observed proton 
and neutron kinetic energy spectra using a simple model 
that included the effects of final-state interactions (FSI) . 
In the model, two types of decay events were considered: 
A7V decay with no FSI (referred to here as no-scatter) 
and decay events where one of the decay nucleons scat- 
ters from one of the spectator nucleons (rescatter). The 
expected no-scatter and rescatter nucleon kinetic energy 
distributions (shown in Fig.[8| were derived from a Monte 
Carlo simulation of the AN interaction. (Also shown is 
the distribution for nucleons from ANN decay, which will 
be discussed later in this section.) For the Monte Carlo 
events, the momenta of the four initial-state baryons were 
generated from a Gaussian distribution with an RMS 
width of 125 MeV/c for the nucleons and 50 MeV/c for 
the A. The distributions were correlated such that the 
sum of the four momenta added to zero, and the energy 
of each spectator nucleon was determined by setting its 
mass to be on-shell. For no-scatter events, the energy 
distribution for the primary decay nucleons was deter- 
mined by choosing a direction for the participating nu- 
cleon isotropically in the A-nucleon center-of-mass frame 
and boosting to the laboratory system. The re-scatter 
events were simulated by choosing a direction for one of 
the nucleons isotropically in the decay-spectator nucleon 
system (s-wave) and boosting to the laboratory frame. 

As in the n~ case above, the expected kinetic energy 
distributions were fit to the observed spectra, and the 
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ity per spectator nucleon, e, and the nucleon-stimulated 
branching fractions, B p and B n , then summing over no- 
scatter and rescatter nucleons yields: 



40 60 80 100 120 140 

Kinetic Energy (MeV) 

FIG. 8: Shown here are the predicted no-scatter and rescatter 
kinetic energy distributions for nucleons from AiV decay (solid 
and dashed lines, respectively) along with the distribution for 
the primary nucleons from ANN decay (dotted line). 



results were combined with the model to extract the 
nucleon-stimulated decay rates. To include the effects of 
the resolution and acceptance of the detector system, the 
expected no-scatter and rescatter distributions of Fig. [8] 
were used as input to the Monte Carlo simulation dis- 
cussed in Sec. |Vl|and then reconstructed using the anal- 
ysis code. The resulting spectra were designated h°(E) 
for the nucleons that did not undergo FSI and h}(E) 
for the rescattered nucleons. (These distributions differ 
for protons and neutrons due to differing energy resolu- 
tions and detection acceptance.) The distributions are in 
units of counts/bin per event and are normalized to the 
detector acceptance. The fits were then performed using 
functions of the form: 



N a {E) = M° a h°a(E) + MX(E), 



no-scatter 



rescatter 



(18) 



where a = p,n and the fit parameters M° and M\ repre- 
sent the acceptance-corrected numbers of no-scatter and 
rescatter nucleons, respectively. The fit parameters M° 
and were related via our model as described below. 

The effects of final-state interactions were incorporated 
by defining a rescatter probability, e, as the probability 
for one of the primary-decay nucleons from AiV decay to 
scatter off of one of the spectator nucleons. The value of 
e was determined from the data as described below. The 
scattering processes were assumed to occur with equal 
probability (i.e., e has the same value for all proton- 
neutron rescatter combinations) . Since the effects of the 
FSI correction were seen to be small, this assumption 
was considered adequate. Assigning a probability to each 
possible final state in terms of the rescatter probabil- 



M° = (l-2e)B p Ngt, (19) 

M° = (1 - 2e)B p N^ + 2(1 - 2e)B n N^, (20) 

Ml = 4e(B p + B n )N^, (21) 

A£ = 4e(B p + B n )Ngl (22) 



Using these expressions, the parameters M^ and M\ were 
eliminated from the fit functions: 

N p (E) = M y p{E ) + ^ € {M a p +M a n )hl{E),{2-i) 

N n (E) = M°hl{E) + ^ e {M° p +Ml)h 1 n {E){2A) 

The three parameters M°, M°, and e were then deter- 
mined by a fit to the data (shown in Fig. |9j. Since the 
fit shapes h°(E) and h}(E) include the effects of the de- 
tector efficiency and resolution, the observed nucleon ki- 
netic energy spectra used for the fits were not efficiency 
corrected. The fits were carried out simultaneously by 
varying the quantity e to obtain a minimum combined 
Xi- Using the fit values, the nucleon-stimulated branch- 
ing fractions were then calculated using Eqs. [19] and [20] 
and the results are listed in Table |IV| The upper lim- 
its for the neutron-stimulated branching fraction and the 
ratio r„/r p were determined by the prescription of Feld- 
man and Cousins [12] , which defines the confidence levels 
in a consistent way while imposing the physical bound- 
ary at zero. The upper limits are given at the 95% CL. 
Also listed are the results for the case of no FSI found by 
fixing the rescatter probability e to zero. 

The model was also extended to include possible con- 
tributions from three-body ANN decays. The expected 
kinetic energy distribution for these nucleons was gen- 
erated using the same Monte Carlo code as for the 
AiV interaction described above, and the resulting dis- 
tribution is shown in Fig. [8] In this case, the mo- 
mentum of the single on-shell spectator was initially 
selected, fixing the momentum and energy of the par- 
ticipating ANN system. Final-state momenta for these 
three baryons were then selected in the system's center- 
of-mass frame according to a phase-space distribution 
generated by the CERNLIB GENBOD code. The nu- 
cleon momenta were then boosted to the laboratory 
frame. 

As in the case of FSI above, the contribution of 
ANN decay, r m b/r tot , was determined by a fit to the 
observed kinetic energy spectra of Fig. [9] Because the 
ANN decay nucleons and rescattered AiV nucleons are 
indistinguishable at the level of statistics of the current 
measurement, the fits were performed with an h}(E) that 
consisted of a linear combination of the ANN and rescat- 
ter distributions. The results for the nucleon-stimulated 
branching fractions are shown in Table [TV] for the cases of 
(i) ANN decay with no FSI for the AiV decays (i.e., e = 
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TABLE IV: Results of the nonmesonic decay rate extraction considering FSI and ANN decay contributions. The decay rates 
are given in units of Ftot, and quoted errors are statistical only. The upper limits given for T n /r to t and F n /F p (listed directly 
beneath the associated quantities) are at the 95% confidence level and include systematic errors. 



No FSI or ANN 



FSI 



ANN 



FSI + ANN 



Systematic Error 



X 2 /d.o.f. 
e 

Fmb/Ftot 
Fp/Ttot 

r„/rtot 

/T tot 



0.78 



0.158 ±0.011 
-0.0066 ± 0.0060 

< 0.018 
0.1510 ±0.0099 
-0.042 ± 0.037 

< 0.10 



0.72 
0.045 ±0.035 

0.160 ±0.015 
-0.0089 ± 0.0075 

< 0.017 
0.151 ±0.013 
-0-056±g : gg 

< 0.10 



0.69 

0.042 ±0.028 
0.134 ±0.026 

_n m qq+0.0089 

u.uioo_ 0071 

< 0.015 
0.162 ±0.020 

-0.100 ±0.068 

< 0.097 



0.71 

O.Q30t°;°2 

0.0171°;°^ 

0.150 ±0.016 

-0 0112+ 0084 
u.uiiz_ 0077 

< 0.017 
0.156 ±0.017 

-0.075 ±0.055 

< 0.098 



±0.0089 
±0.0104 

±0.0096 
±0.063 
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FIG. 9: Shown here are the fits of the expected no-scatter 
and rescatter distributions to the observed (a) proton and (b) 
neutron kinetic energy spectra. The no-scatter contribution 
is indicated by a dashed line, while the rescatter contribution 
is shown as a dotted line. 



0) and (ii) a mixture of ANN decay and FSI with equal 
strengths (i.e., the quantity consisted of a 50/50 
mixture of ANN decay events and rescattered AN decay 
events). From the case of FSI, the upper limit for the con- 
tribution of final-state interactions to the total \He decay 



- A 

rate was found as T^/Ttot = 4e(r mrt /T to t) < 0.11 
(95% CL). Also, the multibaryon decay branching frac- 
tion found for the case of ANN decay with no FSI gives 
an upper limit of T mb /T tot < 0.097 (95% CL) for the 
contribution of the ANN decay process. 



D. Nonmesonic rates from multiple nucleon 
coincidence data 



The AA" interaction produces two energetic nuclcons 
that may be detected in the out-of-beam detector ar- 
rays resulting in both proton neutron (pn) and neutron- 
neutron (nn) coincidence events. The observation of such 
events gives another window onto the nucleon-stimulated 
decay process. In the present experiment, a total of 87 
pn coincidence events and 19 nn coincidence events were 
observed. These events were subject to the same in-beam 
tracking cuts as the excitation spectra of Sec. |IV| along 
with a ±2cr excitation energy cut. The protons and neu- 
trons were identified in the out-of-beam detector arrays 
as described in Sec. IIIII 

The cosine of the separation angle for the two final- 
state nucleons, cosOnn, is shown in Figures [Toj^a) and 
10 'b) for pn and nn coincidence events, respectively. If 
the spectator nucleons are ignored, the two nucleons re- 
sulting from the AA~ interaction would be expected to 
emerge roughly back-to-back with a separation angle 
near 180° (from simple momentum conservation). The 
observed distribution for pn events is strongly peaked 
near the value of 180° with only about 10% of the pn 
events having a separation angle less than 140°. The 
latter events may be indicative of FSI or ANN decay 
contributions. The number of counts is low in the nn 
spectrum, but the back-to-back peak does not seem to 
be as pronounced for these events. The distributions 
for the sum of the kinetic energies of the coincident nu- 
cleons are shown in Fig. [TT] For the AA" interaction, 
the total energy available to the two decay nucleons is 
~166 MeV (i.e., the difference in the A and n masses 
less the binding energy of the ^He hypernucleus). The 
pn kinetic energy distribution of Fig. 11 'a) peaks near 
166 MeV as expected. The observed spread is largely 
consistent with a spread due to the Fermi momentum 
of the baryons within the nucleus coupled with detector 
resolution, and the apparent tail out to lower energies 
may be the result of FSI or ANN decays. (The observed 
spectrum is in qualitative agreement with the results of 
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FIG. 10: The cosine of the nucleon-nucleon separation angle 
is shown for (a) pn and (b) nn coincidence events. Overlaid 
on these plots is the relative geometric acceptance for NN 
coincidence events normalized to the number of counts in each 
spectrum. 



a Monte Carlo simulation which includes such effects.) 
The nn kinetic energy distribution of Fig. |TlJb) , on the 
other hand, does not behave as expected and may be in- 
dicative of nn coincidence events originating from sources 
other than An decays (i.e., some mixture of background 
neutrons and neutrons from Ap or ANN decays). 

The nucleon-stimulatcd branching fractions were ex- 
tracted from the observed pn and nn coincidence events 
by relating the number of each type of coincidence event 
to the branching fractions as: 



N„ 



tot 
HNi 



(25) 



Nun = (P&P£ n + P£)B n g 2 Nt&+N% i , (26) 



where g 2 is the fraction of a Gaussian distribution con- 
tained within ±2er of the mean, Pall is the probability of 
detecting exactly one accidental neutron per event, Pp P 
is the probability of detecting one proton from Ap decay, 
PpP is the probability of detecting one proton and one 
neutron from Ap decay, and P Arl and P^n are the proba- 
bilities of detecting one and two neutrons from An decay, 
respectively. With the exception of Pacl, which was ex- 
tracted from the data, the above detection probabilities 
had to be estimated using the Monte Carlo simulation de- 
scribed in the previous section. Final-state interactions 
were included in the Monte Carlo at the level determined 
for the case of FSI only. 



FIG. 11: The sum of the kinetic energy for the two observed 
nucleons is shown for (a) pn and (b) nn events. The line is 
drawn at the expected value of 166 MeV. 



The term N%> n appearing in Eq. 26 is the expected 



number of nn background events due to accidental neu- 
trons, neutrons from 7r _ absorption, and neutrons from 
proton-stimulated decays. As with the accidental de- 
tection probability above, the detection probabilities for 
neutrons from 7r _ interactions were determined from the 
data with no presumptions about the extent of FSI or 
ANN decays, while the detection probabilities for neu- 
trons from Ap decay were again estimated using the 
Monte Carlo with FSI included. The number of ex- 
pected background nn coincidence events was found to 
be N^ 9 n — 18.3 ±6.1 with the largest contribution arising 
from neutrons from 7r~ absorption (~13 events). Such 
a large background is consistent with the observed nn 
energy-sum spectrum of Fig. |ll[ b). 

The resulting values for the branching fractions are 
listed in Table [V] along with the results from the analysis 
of the single-particle spectra of the previous section for 
the case of FSI only. Also shown is the weighted average 
of the two determinations. 

In addition to the pn and nn events, a handful of pp 
coincidence events were also observed and were used to 
set an upper limit on the App branching fraction. After 
the in-beam tracking cuts and ±2<r excitation energy cut, 
a total of 8 pp coincidence events remained. Although 
these events may be indicative of a App three-body decay 
mode, they could also be Ap decay events for which the 
spectator proton has been detected. As in the pn and 
nn cases, the number of pp events of each type may be 
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TABLE V: Nonmesonic decay rates as determined from the 
single-particle kinetic energy spectra and multiple nucleon co- 
incidence data. The errors given consist of statistical plus 
systematic errors, and upper limits are given at the 95% CL. 
Also shown is the weighted average of the two determinations. 





Single-Particle 


Multiple Coincidence 


Average 


r P /r tot 


0.160 ± 0.018 


0.165 ± 0.028 


0.161 ± 0.015 


r„/r tot 


-0.0089 ± 0.0135 


0.0031 ± 0.0338 


-0.0073 ± 0.0125 




< 0.018 


< 0.070 


< 0.018 


r„ m /r tot 


0.151 ± 0.017 


0.168 ± 0.042 


0.153 ± 0.016 


r„/r p 


-0.056 ± 0.082 


0.019 ± 0.205 


-0.046 ± 0.076 




< 0.11 


< 0.42 


< 0.11 



TABLE VI: The A He decay rates as determined by the 
present analysis. The errors quoted include contributions 
from systematic error, and upper limits are given at the 95% 
CL. 



Mesonic rates 



Nonmesonic rates 



iWr tot 0.564 ± 0.036 iyr tot 
v^/Vtot 0.270 ±0.024 r„/r tot 



r T° +r , r 

r w »/r« 



0.835 ±0.021 
2.09 ±0.29 



i/Ttot 



1*71 /Tp 



0.169 ±0.019 
-0.0032 ±0.0183 

< 0.032 
0.165 ±0.021 

-0.019 ±0.108 

< 0.19 



expressed in terms of the A He branching fractions as: 



function of the form: 



Nif = g 2 s n P± p B n N^ 



pp 

lyApp _ 
pp 



"HNi 



q ± pp -"p- 
2 e n P± pp B, 



(27) 
(28) 



where e q = 0.986 is the measured efficiency of the charged 
trigger, and Pp P pp are the probabilities of detecting 
two protons from Ap and App decay, respectively, and 
B pp is the App branching fraction. The probability of 
detecting two protons from Ap decay is negligibly small 
when interactions with the spectator nucleons are ig- 
nored, but when FSI is included at the level determined 
in Sec. |VII C| it becomes significant. The probabilities 
Pp P and Ppp PP were estimated using the Monte Carlo 



simulation discussed in Sec. VII C including FSI for the 
Ap case. The expected number of back ground pp events, 
N£ p , was then calculated from Eq.|27]as 6.40 ±0.69. As- 
suming Poisson statistics and using the prescription of 
Feldman and Cousins |13j . the upper limit for the num- 
ber of observed pp events originating from App decay was 
found as N^ pp < 9.29 (95% CL for a total number of 8 ob- 
served counts and an expected background of 6 counts) . 
Using Eq. 28 the upper limit for the App branching frac- 



tion was then found as B pp < 0.026. 

As mentioned above, the contribution from Ap decay 
events becomes negligible when final-state interactions 
are ignored. In this case, all 8 observed pp events 
are considered to have originated from App decay, and 
the App branching fraction is found to lie in the range 
0.008 < Bp P < 0.041 at the 95% CL. 



E. Combined \He decay rates 

The final values for the A He branching fractions were 
determined by considering a set of parameters (B^-, 
B^o, Bp, B n , N) representing the 7r~, 7r°, proton- and 
neutron-stimulated branching fractions along with the 
total number of ^He hypernuclei. The optimal values 
for this set of parameters were found by minimizing a \ 2 



2 (Af ff - -B„-NY 
X a 2 _ 



(M^-B^NY 



(Mp-BpN) 2 (M„ - (Bp + 2B n )N)" 



+ 



(M pn -_BpN) 2 , (M, 

T 2 



nil -E$7i N ) 
2 



(29) 



a 



i in 



N 



where the values M a are the acceptance-corrected num- 
bers of event type a with associated error a a . Addition- 
ally, the parameters representing the branching fractions 
were subject to the constraint: 



B„ 



B n + B^~ + B^o — 1, 



(30) 



and any possible contributions from three-body 
ANN decay, which were previously seen to be small, 
were neglected. The minimization was carried out using 
the ROOT MINUIT minimization routine, and the 



results are listed in Table VI (in units of T to t) 



VIII. DISCUSSION 



The ^He decay rates determined in the present work 



are listed in Table VII (in units of Ta). Also listed are the 
results of another measurement of ^He decay performed 
by Outa et al [2] along with the results of some earlier 
bubble chamber experiments [5J [T31 [TS] ■ The values taken 
from reference [8] are the result of a reanalysis of various 
older results performed by Schumacher. The decay rates 
found in this work are seen to be in good agreement with 
past measurements. 

The results reported for the experiment of Outa et al 
were published in 1998, some 8 years after the running 
of our experiment. Their experiment created ^He by 
the 4 He(K~ cd ,7r~) reaction and detected the 7r°, ir~ , 
and protons from the subsequent weak decay. Unlike 



14 



our experiment, their apparatus was not able to de- 
tect neutrons, requiring them to determine the neutron- 
stimulated rate by subtraction. Our results for the 
neutron-stimulated rate and neutron to proton ratio, 
r n /rA < 0.035 and r„/r p < 0.19, provide stricter upper 
limits as compared to T„/rA < 0.09 and T n /T p < 0.60 
for Outa et al (all at the 95% CL). The upper limits 
for Outa et al given above were derived from the results 



listed in Table VII using the same method as for our re- 
sults. There is also a slight difference in the mesonic de- 
cay rates, where our T^o /T n - ratio is seen to be in better 
agreement with the earlier bubble chamber results. 

Selected theoretical investigations of the 
^He nonmcsonic partial decay rates are shown in 
Table VIII These OPE (one pion exchange) models 



incorporate a quark-based model for short range in- 
teractions and are currently considered to be the most 
successful of the OPE-based models. The results of 
Jun [T7] and Inoue et al [5] are shown for OPE only and 
OPE + quark interaction mechanism. The work of Jun 
utilizes a phenomenological 4-baryon point interaction 
model to describe the short range quark interactions, 
while Inoue et al uses a direct quark (DQ) interaction 
based on the effective weak Hamiltonian of Gilman and 
Wise [18]. The calculation of Sasaki et al [6] employs 
the direct quark mechanism of Inoue et al coupled with 
a 7r + K ± a exchange potential. Unlike previous models, 
this model is able to reproduce the value of the proton 
asymmetry observed in polarized ^Hc nonmesonic 
decay. All of the above models adequately reproduce 
the observed ^He nonmesonic decay rates. 

In regards to the A I = | rule, it is interesting to note 
that the direct quark models predict a significant con- 
tribution from Al = | transitions arising from the ef- 
fective weak Hamiltonian used, while models based on 
the phenomenological 4-baryon point interaction are able 
to reproduce the experimental results with little or no 
Al = | component. All of these models enforce the 
Al = 1 rule for the meson-exchange component. 



TABLE VII: Measured properties of the A He hypernucleus. 
Listed here are the results of the present analysis along with 
the results of another recent experiment on A He decay [14) 
and some older results from bubble chamber experiments [101 
1151 116] . The errors quoted for this analysis include contribu- 
tions from systematic error. 



This work 



Outa et al [14: Earlier results 



Ttot /Ta 

rv/r A 
iV /r A 
ivAV 

r.+AV 

r p /r A 

r n /r A 

/Ta 

T nm /r^Tj- — 



1.07 ±0.11 1.03±g io 

0.604 ±0.073 0.53 ±0.07 

0.289 ±0.039 0.33 ±0.05 

2.09 ± 0.29 1.59 ±0.20 



0.180 ±0.028 

< 0.035 
0.177 ±0.029 
0.61 ±0.10 

< 0.19 



0.16 ±0.02 

o m+ 004 

u - ul -o.oi 
0.17 ±0.05 

0.51 ±0.16 

06+ 0,28 
u - uo -0.06 



1.15 ±0.48 8 



2.49 ±0.34 [15] 
2.20 ±0.39 \W\ 
0.043 ±0.017 HO] 



0.56 ±0.09 [IS] 
0.40 ±0.15 [15] 
0.29 ±0.13 QS] 



TABLE VIII: Results of selected theoretical investigations of 
the iHe partial decay rates. The partial rates listed here are 
given in units of the free A decay width Ta. 





Jun [I7] 


Inoue [5] 


Sasaki [6] 




OPE 


OPE±4BPI 


OPE 


OPE±DQ 


DQ± 


r„/r A 


0.0008 


0.0373 


0.009 


0.038 


0.015 


r p /r A 


0.1478 


0.1602 


0.145 


0.214 


0.169 


Trim /Ta 


0.1486 


0.1975 


0.154 


0.253 


0.184 




0.005 


0.23 


0.061 


0.178 


0.091 



A. FSI, three-body ANN decays, and interference 
effects 

Using the simple model for nonmesonic decay intro- 
duced in Sec. | VII C] the contributions of final-state inter- 
actions and three-body ANN decays were studied utiliz- 
ing the single-particle proton and neutron kinetic energy 
spectra, and upper limits for each case were determined 
as: 

< 0.11, and — < 0.097, (31) 
Ttot Ttot 

for final-state interactions and ANN decays, respectively 
(95% CL). Within the model used here, the effects of FSI 
and ANN decays on pertinent quantities were seen to be 
small and, thus, have little effect on the final interpreta- 
tion of the data. 



In the analysis of the nonmesonic decay rates, the var- 
ious decay modes have been assumed to add incoher- 
ently, neglecting possible quantum interference effects. 
It has been suggested [IHJ [20] that such interference ef- 
fects may be the origin of the so-called r„/r p puzzle, 
and Garbarino et al |2Q] point out that any interference 
effects would be less pronounced for multiple coincidence 
measurements, providing a cleaner extraction of the non- 
mesonic decay rates. For the case of ^He presented here, 
the results of the two determinations (listed in Table [v]) 
are consistent , suggesting that interference effects are not 
significant at this level of statistics. A detailed study of 
the angular correlations in the multiple nucleon coinci- 
dence data might shed more light on the contributions of 
interference effects as well as FSI and ANN decays, but 
the number of coincidence events in the present data set 
is too limited to provide anything conclusive. 
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B. Phenomenological analysis of s-shell hypernuclei 

As mentioned in the introduction, measurements of 
the nonmcsonic decay rates for the s-shell hypernucleus 
A H are currently quite limited. The most recent results 
from Outa et al jTJj give a value for the total nonmesonic 
rate of r nm /r A ( A H) = 0.17 ±0.11, but no measurements 
of the proton- and neutron-stimulated partial rates have 
been made. Using the present results for A He, these par- 
tial rates can be calculated in the context of the pre- 
viously described phenomenological model of Block and 
Dalitz. 

From Eq. [7] the upper limit of r„/r A ( A He) < 0.035 
determined for the current measurement implies an up- 
per limit of r p /r A ( A H) < 0.017 for the A H proton- 
stimulated rate under the assumption of pure AI = i 
transitions. A measurement of the A H proton-stimulated 
rate greater than this value would indicate a violation 
of the AI = \ rule for the AN weak interaction. The 
A H neutron-stimulated decay rate was found from Eq. [8] 
as r„/r A ( A H) = 0.081 ± 0.023. This calculation utilizes 
the result of Outa et al for the ratio of the neutron- to 
proton-stimulated rates for A He, r„/r p ( A He) = 0.45 ± 
0.11 ± 0.03 [9], but makes no assumptions regarding the 
AI = ^ rule. Using this result, an alternate determina- 
tion of the A H proton-stimulated rate can be made from 
the result of Outa et al for the total nonmesonic decay 
width of A H: 

r p ( A H) - r„ m ( A H) - r n ( A H) = (0.09 ± o.n)r A , (32) 



which, again, is independent of any assumptions about 
the isospin nature of the AN interaction. Because of 
the large error, this result is consistent with the upper 
limit found above. Thus, at the present level of statistics 
for the A H measurement, nothing conclusive can be said 
about the AI = | rule for the AN weak interaction. 

The values of the A H nonmesonic decay rates derived 
from the above discussion represent the most current 
knowledge (although indirect) of the A H proton- and 
neutron-stimulated partial rates. The upper limit deter- 
mined for the A H proton-stimulated decay rate suggests 
a direct measurement of this decay mode as an easy test 
of the AI = I rule for nucleon-stimulated decays (i.e., a 
measurement in excess of this upper limit would show a 
clear violation of the AI = g rule for the AN weak inter- 
action) . Measurements of the A H nonmesonic decay rates 
along with more precise determinations of the A He and 
A Hc nonmesonic decay rates would help to further deter- 
mine the isospin structure of the AiV weak interaction. 
Multiple nucleon coincidence studies can potentially pro- 
vide the cleanest extraction of the nucleon-stimulated 
partial rates and would be a worthwhile avenue of future 
research. 
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